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ABSTRACT 

We compare three analytical prescriptions for merger times available from the lit- 
erature to simulations of isolated mergers. We probe three different redshifts, and 
several halo concentrations, mass ratios, orbital circularities and orbital energies of 
the satellite. We find that prescriptions available in the literature significantly under- 
predict long timescales for mergers at high redshift. We argue that these results have 
not been highlighted previously either because the evolution of halo concentration 
of satellite galaxies has been neglected (in previous isolated merger simulations), or 
because long merger times and mergers with high initial orbital circularities are under- 
represented (for prescriptions based on cosmological simulations). Motivated by the 
evolution of halo concentration at fixed mass, an explicit dependence on redshift added 



mod 
merger 



(z) = (1 + z)°-^H- 



to the prescription based on isolated mergers gives a 



as t 

significant improvement in the predicted merger times up to ~20 tdyn in the redshift 
range 0^z^2. When this modified prescription is used to compute galaxy stellar mass 
functions, we find that it leads up to a 25 per cent increase in the number of low mass 
galaxies surviving at z=0, and a 10 per cent increase for more massive galaxies. This 
worsen the known over-prediction in the number of low mass galaxies by hierarchical 
models of galaxy formation. 

Key words: galaxies: evolution - galaxies: structure - galaxies: kinematics and 
dynamics - galaxies: interactions - methods: N-body simulations 



1 INTRODUCTION 

Galaxies orbiting in dense environments, such as groups or 
clusters, suffer a continuous loss of energy and angular mo- 
mentum under the effect of dynamical friction against the 
medium of their parent halo. Because of this, galaxies on 
bound orbits spiral-in towards the densest regions of their lo- 
cal environment on a given timescale. This merger timescale 
tmerger is usually modelled as a function of the dynamical 
timescale of the main halo, fdyn, the galaxy- main halo mass 
ratio, and the orbital energy and circularity of the infalling 
galaxy. All of these quantities are considered at the time a 
galaxy crosses the virial radius of the main halo along its in- 
falling orbit. Merger timescales are defined as the time until 
the galaxy has lost a significant fraction of either its ini- 
tial mass via tidal stripping, or of its initial orbital angular 
momentum. 

Accurate estimates of how long a galaxy survives within 
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a halo while being affected by dynamical friction are funda- 
mental for theoretical studies of galaxy evolution. Indeed, 
merger times play a key role in the evolution of galaxies' 
stel lar m asses, mor phologies, colours , and gas content (e.g. 
iCox et al.. 2008: De Lucia et al.ll201(]| . [20111 . 

In recent years, iBovlan-Kolchin et all jo doj, [Boj) 
iJiang et all l|2008l . [joj ). and iMcCavana eTall (|2012l . llVld ) 

have provided different prescriptions to estimate merger 
timescales. All three prescri ptions are variati o ns of the ana- 
lytic description derived by I Chandrasekhaij (| 19431 ') for the 
drag force suffered by a point mass object as it moves 
through a uniform background mediu m of less massive par- 
ticles (see iBinnev fc Tremainel Il987l ). The aforementioned 
prescriptions are obtained either by simulating single iso- 
lated mergers at z=0 that probe a given parameter space, 
or by collecting a sample of mergers from cosmological simu- 
lations within a given redshift range. In this paper, we com- 
pare these prescriptions with controlled simulations of iso- 
lated galaxies being accreted onto a group-like halo at three 
redshift epochs. We explore several halo concentrations in 
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Table 1. Properties of group environments and galaxies. 





"z = 0" 


"z = l" 


"z=2" 




Group DM halo 










Virial mass 


9.9 


9.9 


9.9 


(xIO^^Mq) 


Virial radius 


555.94 


329.19 


226.09 


(kpc) 


Concentration 


9.74 


4.87 


3.25 




Circular velocity 


276.97 


360.07 


434.62 


(km s^^) 


Number particles 


5.5 


5.5 


5.5 


{X10=) 


Softening 


0.55 


0.32 


0.22 


(kpc) 


Group stellar spheroid 








(xl0"AfQ) 


Mass 


1 


1 


1 


Scale radius 


3.24 


1.91 


1.31 


(kpc) 


Number particles 


2.5 


2.5 


2.5 


(xlO^) 


Softening 


0.1 


0.06 


0.04 


(kpc) 


Galaxy DM halo 








(xlO^AfQ) 


Virial mass 


2.97-39.6 


2.47-34.6 


2.97-9.9 


Virial radius 


172.7-409.6 


96.2-231.9 


70.2-104.9 


(kpc) 


Concentration 


15.36-10.97 


7.87-5.58 


5.12-4.38 




Number particles 


2.5 


2.5 


2.5 


(XlO^) 


Softening 


0.35 


0.26 


0.19 


(kpc) 


Galaxy stellar disc 










Mass 


2.8 


1.42 


0.72 


(XlO^^Afg) 


Scale-length 


3.5 


1.65 


0.9 


(kpc) 


Number particles 


5 


5 


5 


(XlO*) 


Softening 


0.05 


0.012 


0.007 


(kpc) 



both systems, merger mass ratios, and orbital parameters 
of accreted galaxies. Our main goal is to determine whether 
the implicit dependency of these prescriptions on redshift 
(via tdyn) is enough to account for relevant properties of 
haloes that are known to evolve with cosmic time, such as 
halo concentration. 

The layout of this paper is as follows: Section [2] de- 
scribes the set-up of our experiments; Section [3] describes 
our results, comparing them to models from the literature; 
in Section |4] we discuss our results and in Section [5] we give 
our conclusions. 



2 SET-UP OF NUMERICAL EXPERIMENTS 

We have carried out 50 simulations of isolated mergers be- 
tween a single galaxy and a larger main hal o, in order to 
quantify their merger timescales. Similarly to IBOSI . our ba- 
sic strategy is to release a single disc galaxy at a time, on 
a bound orbit, from the virial radius of the main halo and 
study the evolution of its mass content and orbital angu- 
lar momentum as it is dragged towards the centre of the 
halo. Our simulations cover the time span since a galaxy is 
released from the virial radius of the main halo until it is ei- 
ther disrupted or it has exhausted its initial orbital angular 
momentum (see Section r3.1[) . 

Each main halo is modelled as a A^'-body self-consistent 
multi-component system, where the DM halo follows a NFW 
density profile (Navarro et al. 1997i ): 

Phalo(7-) = , , -1—1-^, (1) 

where ps is a characteristic scale density and a scale ra- 
dius. Additionally, a spherical stellar component, resembling 
a central galaxy, is located at the halo centre, and its mass 
follows a Ifernquist density profile (|Hernauist|[T99Q ) : 



where M* is the stellar mass and a, is the scale radius. Sim- 
ilarly, each galaxy is also modelled as a self-consistent multi- 
component system, formed by a stellar disc embedded in a 
DM halo. The stellar mass in discs follows an exponential 
density profile: 

Pdisc(fl,2) = J^^^" exp ( --^ ) sech^ ' 

where A/disc is the disc mass, Rq is the exponential scale- 
length, and 2D is the exponential scale-height. The DM halo 
of galaxies is also assumed to follow a NFW profile. All DM 
haloes in our simulations are initially spherical, do not ro- 
tate, and the structure of their inner region has been adia- 
batically contracted to account for the growth of the bary- 
onic component (see IVillalobos fc HelmilboOSl . for details). 

The initial orbital parameters of discs are chosen to 
be consistent with distributions of orbital parameters of 
infalling su bstructures, o b tained from c osmological simula- 
tions (e.g.. iBensonI l2005l : IWetzell l 201ll ). In this paper, we 
focus on the most likely orbital circularity of infalling sub- 
structures, and on the extreme values of the distributions. 

The initial conditions of our isolated mergers are placed 
in a simplified context of three different redshift epochs, 
"2=0", "2=1" and "2=2". In this context, the initial struc- 
ture of DM haloes is defined by both the virial over-density 
Avir(2) (iBrvan fc Normanlll99^ and the halo concentration 
c(Mvir, 2)0 Note that during the simulations, a main halo 
only interacts with a single disc galaxy, and does not accrete 
additional mass. As for the evolu tion of halo concent ration, 
we adopt c(Mvir, 2) oc (1-1- 2)"^ (|Bullock et al.ll200ll '). 

At all three redshifts, main haloes have a mass of 
W^^M(^). Thi s choic e is b ased on the mass range reported by 
iMcGee et al] (|2009l ') and iDe Lucia et al.l (|2012l '). where sig- 
nificant e nvironmental effects on galaxies must take place 
(see also iBerrier et al.l :2009). We have kept fixed both 
the mass and scale-lengths of stellar discs, while covering 
a range of masses and radii of the DM haloes in which 
they are embedded (see below). The stellar-to-DM mass ra- 
tios of disc galaxies are consistent with stellar-to-DM mass 
relations from both observations and theoretical studie s 
(jBehroozi et al.ll2010l : IGuo et aLllioiol : iMoster et aDl2010l ). 
Table [T] lists the structural parameters of main haloes and 
disc galaxies. Table [5] provides the range of the parameter 
space covered by our experiments. These parameters delimit 
the validity of our results. Note that our experiments do not 
include cases where Chost ~ Csat which would be relevant in 
the case of major mergers. 

We have carried out test simulations to ensure that our 
results are not affected by numerical resolution. All simula- 
tions were run using GADGET-3 (tSpringel,2005i '). 



3 RESULTS 

3.1 Definition of "merger" 

As a galaxy orbits within a halo, it gradually loses orbital 
energy and orbital angular momentum due to dynamical 
friction (mainly) against the medium of the larger halo. As 
its orbit decays towards the densest region of the halo, a 

1 We assume Qm,0 = 0.3, Ha = 0.7 and Hq = 70 km 8"^ Mpc~l. 
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Table 2. Ranges of parameters explored by the experiments. 



z 


-MhoBt/Msat 




rc(B)/rvir 






tatara 


(1) 


(2) 


(3) 


(4) 


(5) 


(B) 


(7) 





2.5-33.33 


0.21-0.91 


0.91-1.94 


4.8-26.5 


4.8-28.5 


5.3-29 


1 


2.86-40 


0.21-0.9 


1.01-1.97 


3.2-22 


3.7-26 


4.2-26 


2 


10-33.33 


0.55 


1.06-1.16 


2.5-10 


3.5-22 


4-24 



(1) Redshift. (2) Mass ratio of each merger. (3) Initial orbital circularity of 
the satellite, in units of ^c{^vir)- (4): Initial orbital energy of the satellite, 
expressed in terms of the radius of a circular orbit with the same orbital 
energy £/, in units of 7Vir. (5): merger time computed when only 5 per 
cent of the DM halo of a satellite remains bound, in Gyr. (6): merger time 
computed when only 5 per cent of the initial orbital angular momentum of 
a satellite is retained, in Gyr. (7): merger time computed when only 5 per 
cent of the stellar content of a satellite remains bound, in Gyr. 

galaxy might also undergo mass loss via tidal stripping, espe- 
cially at pericentric passages where tidal forces are stronger. 
Previous studies have used several criteria to define when a 
galajcy has "merged" with a larger halo. A galaxy is consid- 
ered merged if either: (i) its remaining bound DM mass is 
5 per cent of its initial value; or, (ii) the remaining orbital 
angular momentum is 5 per cent of its initial value; or, (iii) 
the remaining bound stellar mass is 5 per cent of its initial 
value. In general, all three definitions yield similar merger 
times, where usually tum <tj< istars- In the rest of the pa- 
per, we adopt the definition of merger time as the minimum 
between tj and fstars- This cor respo nds to tj in most of our 
experiments (e.g. see Fig. 1 in IBOSI ) . 

3.2 Comparison to previous models of merger 
times 

We compare the merger times obtained from our simula- 
tions to th e foll owin g mod els of merger timescales from the 
literature. IBOSI and I Ml 3 . using isolated and cosmological 
simulations respectively, find: 



= A- 



(Af host /Ms, 



(4) 



tdyn ln(l +M"host/Msat) 

with the following best fitting parameters 
(A, B,C,D)=(0 .216.1.3.1.9.1.0) forlBOa and (0.9,1.0,0.6,0.1) 
for iMll iJOSl . using cosmological hydro-dynamical simula- 
tions, find: 

tmerger _ 0.947?°-^ + 0.6 Mhost 1 . . 

tdyn ^ 2X0.43 Msat ln[l-H(Mhost/Msat)]' ^' 

In all three models, tdyn ~ '■vir/V'c(?'vir) = 
y/2/ /S.{z)/ H{z), where the over-density A(2:) and Hubble 
constant H{z) vary as a function of redshift. Note that dy- 
namical timescales do not depend on either halo mass or 
virial radius, but only on the cosmological parameters cho- 
sen. For the parameters used in this paper, the dynami- 
cal timescales in our simulations are fdyn(z=0) = 1.9644 Gyr, 
tdyn(2=l)=0.895 Gyr, and tdyn(2=2)=0.5095 Gyr. 

Figure [T ] (lef t) sho ws the predicted merger times from 
the lBOSl . IJOSI and lMl j presc riptio ns applied to our mergers. 
For experiments at "2=0", iBOa offer accurate predictions 
at least up to ~15 tdyn, independently of the merger defi- 
nition used (Section 13. ip . An analysis of predicted merger 
timescal es as a function of the initial orbit of galaxies shows 
that the IBOSI model tends to slightly overestimate (under- 
estimate) tmerger for galaxi es wi th lower (higher) initial cir- 
cularities (Figure [2|. Also, [bO^ appears to be fine-tuned to 
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Figure 1. Left column: predictions of the [bo3 . Ijoi . ImiI mod- 
els applied to our experiments at "2=0" (black), "2=1" (red), 
and "2=2" (green). Arrows indicate the improved accuracy of 
the prediction for selected experiments at "2=1", when mergers 
are resimulated with a higher initial concentration of the galax- 
ies' haloes. Right column: predictions of the respective modified 



models, 4" 



1.(2) — (1 ~t~ 2)'^tmcrger. 



accurately predict merger timescales for galaxies infa lling 
with the most likely orbital circularities (77~0.6). The IBOSI 
model, applied to our experiments at "2=0", exhibits an 
overall scatter of ±20 per cent across the range of explored 
circularities, decreasing to ±10 per cent over different mass 
ratios at a fixed circularity. The high accuracy of the IBOSI 
prescription is not surprising, since it was obtained from 
simulations of isolated mergers at z=Q, similar to those pre- 
sented in this paper. However, for experiments at "z=l" 
and "z=2", a comparable level of accuracy in the predic- 
tions is only rea ched up to ~8 tdyn (Figure [TJ top-left). Be- 
yond that limit, IBOSI systematically under-predicts merger 
times. Larger differences are observed for experiments with 
long er merger times and for those at higher redshift. As IBOSI . 
IJOSI and |M12| also significantly under-predict longer merger 
times, and increasingly so for mergers at higher redshift 
when compared to our controlled merger simulations. This 
implies that the implicit redshift dependence in all three 
models, from tdyn(-z), is not enough to account for the mea- 
sured evolution of merger timescales. 



4 A. Villalobos, G. De Lucia, S. Weinmann, S. Borgani, and G. Murante 



9 1.5 



0.5 



r» B08 


1 1 1 1 1 1 


1 'a ' ' _ 


- □ JOB 




□ 

A 


: A M12 


A 


□ 

A _ 


: 

... 


A 
° A 
■ □ A 

M-B-^-- 


□ A - 
□ 




B A 

" A 

□ 

□ 

1 < , , 1 , 





0.2 



0.4 0.6 
circularity 17 



0.8 



Figu re 2. Comparison between predictions of the lBOa . IJOSI and 
IMI2I models applied to our experiments at "2=0" , as a function 



of the initial orbital circularity of satellites. Predictions from lBOSi 
and [m12I are shown with offsets r)-0.05 and 7?+0.05, respectively, 
for clarity. 
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Figure 3. Left: difference between galaxy stellar mass functions 
using and models, averaged over 10 haloes in 



merger v .' .^^^^ merger 

each mass range: 5 X IQl^MfTi, IO^^'Mq and 5 X 10^'^ M, 



at 2=0. 



Right: averaged fractional increase (in each galaxy stellar mass 
bin) associated to t'^°f (z), with respect to the tj^orger rnodel. 



In general, we find that the|joi and lM12l models yield 
similar predictions for our experiments at "2=0", although 
showing a significantly larger scatter in compariso n to the 
BOS model (Figure [2}. In spite of the larger scatter, [jOSi and 
Mia do predict accurately the median tmergor (for mergers 
of different mass ratios) for experiments with low and most 
likely initial orbital circularities. This is possibly a conse- 
quence of the fact that, in both studies, mergers are ex- 
tracted from cosmological simulations, and so include a large 
number of radial mergers. Both models, however, under- 
predict the merger timescales in case of high circularities, 
which are associated to long merger times, showin g an off- 
setof '^SS per cent. The scatter associated to both lJOSl and 
|M12| models is ±35 per cent, for mergers with low and most 
likely initial orbital circ ularit i es. 

We argue that the IBOSI . Ijo3 and IM12| models under- 
predict long merger times at high redshift as a consequence 
of both neglecting the evolution of halo concentration in 
satellite, and under-sampling of long merger times. 



namical friction formula (jChandrasekhaill 1943l l . as they are 
applied to mergers of extended objects. 

Motivated by the inverse proportionality between 
^merger and c iu our experiments, we introduce a modification 
to the IBOSI model in the form C?r1;er(-z) = (1 + 2)"ime?ger. 
We find that the inclusion of this modification, with a=0.44, 
offers an improvement in the predicted merger timescales at 
least up to ~20 tdyn for mergers between 0^2^2 (Figure [TJ 
top-right) Note however that a factor ~2 varia tion in con- 
centration at fixed halo mass (jNeto et al.ll2007h . implies a 
~1.4 variati on in tme rger . 

In both|joi and lMll the evolution of halo concentra- 
tion is aheady included in their cosmological simulations by 
construction. No details about the concentration of haloes 
are however discussed in these studies. 



3.2.1 Evolution of halo concentration in satellites 

In order to study the effect of changes in the density profile 
of satellite haloes, we have resimulated three experiments 
at "2=1", increasing the concentration of their DM haloes. 
Specifically, we have given to those DM haloes a higher con- 
centration as haloes of the same mass would have at "2=0" . 
Arro ws in Figure [T] (top-left) show that predictions from the 
IBOSI model become significantly more accurate in the case 
of mergers with more concentrated DM haloes. 

A satellite with a higher concentration is expected to be 
more resilient against tidal disruption. Therefore, it will re- 
tain more of its mass, leading to a shorter merging timescale. 
Since the IBOSI model is obtained from isolated mergers at 
2=0, satellite haloes are assumed to have a too high concen- 
tration when the model is applied to mergers at 2>0, which 
leads to shorter predicted merger times. 

This result strongly suggests that the systernatic under- 
prediction of longer merger timescales by the IBOSI model 
comes as a consequence of the evolution of halo concentra- 
tion in our experiments [c(Mvir,2) oc (1 -I- z)~^ , Section [2], 
which is not included in their simulations. Note that the evo- 
lution of halo concentration also puts a strain on the validity 
of merger times models derived from the Chandrasekhar dy- 



3.2.2 Under-samplmg of long merger times 

Even though cosmological simulations are statistically ro- 
bust and provide a realistic context for studies of merger 
timescales, in general both long timescales (>15 tdyn) 
and mergers with high circularities are severely under- 
represented. For instance, only mergers completed b y 2=0 
and with orbital pericentres ^ rvir are considered in IM12| . 
As a consequence, this likely leads to inaccurate predictions 
of long merger timescales. Even though the motiv atio n be- 
hind our modification does not apply to the IJOSI and IM12| 
models (since evolution of halo concentration is included in 
their models by construction), our proposed modification 
does partially alleviate the disagreement between their pre- 
dictions and our simulations, as shown in Figure[T] Note that 
IBOSI might also be affected by under-sampling given the rel- 
atively low number of long mergers their fitting formula is 
based on. 



^ The time dependence in this modification comes only indirectly 
through the median concentration-redshift relation. 
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4 DISCUSSION 

We study the effect of our modification to the IBOSI model 
on the galaxy stellar mass functions for galaxies residing in 
haloes of different mass. 

The galaxy stellar mass functions are com puted by 
select ing from the Millennium simulation (iSpringel et alj 
120051 ) 30 haloes in three mass ranges (5 x 10^^ Mq, lO" 
Mq, and 5 x 10^* Mq at 2=0). Following their main pro- 
genitors back in time, we extract all information about the 
galaxies residing in each halo and in its main progenitors, 
taking advantage of pub licly available galaxy catalogues 
(|De Lucia fc Blaizotll2007l ). For each galaxy, we then store 
both its mass and that of its host halo at accretion time 
(i.e. the snapshot before it becomes a satellite for the first 
time). In the procedure, we take care of avoiding counting 
the same galaxy more than once. An orbital circularity is as- 
signed randomly to each galaxy f r om d istributions of orbital 
parameters obtained bv iBenso 3 (|2005l ). Then, both t^S,%cr 
and tmorger(2) prescriptions are applied to each galaxy and 
its host. Finally, it is assumed that a galaxy has survived 
by z=0, if imerger>iLB(«accrotion), whcrc ti,B is the look-back 
time at the redshift of accretion. As a first approximation, 
we assume that all surviving galaxies conserve their stellar 
mass since they were accreted. We also assume that the stel- 
lar mass of galaxies that do not survive is added to either 
the diffuse stellar component of the main halo, or to the 
central galaxy. This is clearly a strong simplification, since 
mergers between satellite galaxies within a main halo could 
also take place, altering the intermediate-mass region of the 
mass functions. We plan to study in a future work the con- 
sequences of the proposed modified model in the context of 
a more realistic galaxy formation model. 

Figure [3] (left) shows the difference between mass func- 
tions using tm?J^gor(-z) ^ud imefger modcls, avcragcd over 10 
haloes in each mass range. We find that the modification 
introduced leads to a larger number of lower mass satel- 
lites surviving at z=0, across all explored mass ranges. This 
comes as a natural consequence of the under-prediction of 
long merger timesca les (a ssociated to low mass satellites) at 
high redshift by the lBOa model, as seen in Figure [T] (left). 

Figure [3] (right) shows the fractional increase in the 
number of satellites within each galaxy stellar mass bin, due 
to the use of fm°rger(2) with respect to the tmorgcr model. 
We find that satellites on the low mass end present the 
most significant fractional increase: ~25 per cent consider- 
ing haloes in all mass ranges. On the other hand, more mas- 
sive satellites show a rather constant (within the standard 
errors) fractional increase of ~10 per cent. Previous stud- 
ies comparing observational data to predictions from semi- 
analytic models of galaxy evolution have found that these 
over-predict the num ber of faint galaxies jWeinmann et all 
l2006l : lLiu et"al]|2010l ) . Our results suggest that this problem 
would worsen if more realistic models of merger timescales 
are employed. 

By assuming that all the stellar content of merged 
galaxies contributes to the mass of central galaxies, we find 
that t^tfser{z) leads to a mass reduction of ~10 per cent for 
the central galaxy, with respect to tme?gcr- This is because 
there is less contribution from merged satellites, which would 
have longer survival times. This calculation assumes that the 



initial mass of the central galaxy (i.e., before mergers take 
place) is 1 per cent that of the final halo mass at z=0. 

Similar results are found when we use modified versions 
of both thc iJOS, and .MlZ models instead of iB08 . 



5 SUMMARY AND CONCLUSIONS 

In this work, we compare predictions from three models of 
merger times available in the literature to simulations of iso- 
lated mergers. We generate 50 simulations of single mergers 
between a satellite galaxy and a main halo at three redshift 
epochs, studying the evolution of the galaxy mass content 
and orbital angular momentum, as it is affected by dynam- 
ical friction. We probe a parameter space of different halo 
concentrations, merger mass ratios, orbital circularities, and 
orbital energies of galaxies. 

We find that the implicit dependency on redshift in the 
models is not enough to account for variations as a function 
of redshift in our simulations. In particular, we find that 
prescriptions available in the literature significantly under- 
predict long timescales for mergers at high redshift. In a 
prescription derived from simulations of isolated mergers, 
this is found to be caused mainly by the lack of an explicit 
treatment of the evolution of halo concentration in satellite 
galaxies. On the other hand, in prescriptions derived from 
cosmological simulations, the disagreement is likely due to 
the fact that long merger times, as well as mergers with high 
initial orbital circularities, are under-represented. 

Motivated by the effect of the evolution of halo con- 
centration of satellites, we introduce a modification to the 
model derived from isolated mergers in the form fm°J^gcr (z) = 
(1 4- 2:)"tmcrgor. With a=0.44, the prescription improves sig- 
nificantly the predictions of merger times up to ^--^20 tdyn 
for mergers between 0^z^2. We estimate that our proposed 
modification can lead up to a 25 per cent increase in the 
number of low mass galaxies in massive haloes, and a 10 per 
cent increase in the number of the most massive galaxies. 
This would worsen the disagreement between observations 
and predictions from semi-analytic models of galaxy evolu- 
tion, as found in previous studies. 

Precise predictions of merger timescales are a key in- 
gredient in models of galaxy evolution. In a future work, we 
plan to investigate in detail the influence of our proposed 
modification, in the context of a realistic galaxy formation 
model coupled to A*'-body cosmological simulations. 
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